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ABSTRACT 

Thermonuclear, or Type la supernovae (SNe la), originate from the explosion of carbon-oxygen 
white dwarfs, and serve as standardizable cosmological candles. However, despite their importance, 
the nature of the progenitor systems that give rise to SNe la has not been hitherto elucidated. Obser¬ 
vational evidence favors the double-degenerate channel, in which merging white dwarf binaries lead 
to SNe la. Eurthermore, significant discrepancies exist between observations and theory, and to date, 
there has been no self-consistent merger model that yields a SNe la. Here we show that a spiral mode 
instability in the accretion disk formed during a binary white dwarf merger leads to a detonation on 
a dynamical timescale. This mechanism sheds light on how white dwarf mergers may frequently yield 
SNe la. 

Subject headings: supernovae: general — hydrodynamics — white dwarfs 


1. INTRODUCTION 

It has long been suggested that a white dwarf in a bi¬ 
nary system - with either another white d warf, through 
the so-calle d double-degenerate channel (Iben & Tu- 
tukov||1984|), or a main-sequence or red giaiit companion 
throu gh the single-degenerate channel (Whelan & Iben 
19731 - could give rise to a Type la supernova (SN la) 
event. Yet, despite the recent breakthroug h detections of 
relati vely nearby SNe la, including 2 01 Ife (Nugent et al. 
20111 and 2014J |Zheng et al. 20141, no normal SNe la 
stellar progenitor has ever beeii directly identified. Over¬ 
all, the challenging, unsolved stellar progenitor problem 
for SNe la stands in marked contrast to the case of core¬ 
collapse supernovae, whose stellar progenitors are highly 
luminous, massive stars. 

In light of the challenges to the single-degenerate chan¬ 
nel of SNe la, increasing attention is being paid to the 
double-degenerate channel. In contrast to the single¬ 
degenerate channel, the double-degenerate channel offers 
natural solutions to the stellar progenitor problem, as 
well as to other problem s: the delay-time distribution 
(Maoz & Badenes 20101, the observed absence of dense 
circumstella r gas surrounding the progenitor (Badenes 
et al. 20071, and the absence of the surviving binary 
member star in SNe la remants (Kerzendorf et al. ]M4| ) 

However, the double-degenerate channel m turn 
presents its own set of challenges. Erom an observa¬ 
tional and a theoretical perspective, one of the most out¬ 
standing challenges is to understand the conditions under 
which a merging white dwarf binary system yields a det- 
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onation of the underlying carbon-oxygen white dwarf, 
and hence a SNe la. On this issue, however, there cur¬ 
rently exists a significant discrepancy between observa¬ 
tions and theory. On the one hand, if double-degenerate 
mergers must account for the observed SNe la rate, ob¬ 
servations appear to require that the majority of mergers, 
incl uding sub-Chandrasekh ar mass systems, yield SNe 
la (Badenes & Maoz 20121. In contrast, the most ex- 
tensive numerical studies completed to date suggest that 
essentially no carbon-oxygen white dwarf binary systems 
self-consistently achieve det onation condition s during the 
merg er phase of evolution (Zhu et al. 2013 |Dan et al. 
20141. Simulations of the detonation phase of the so- 
called “violent mergers” of the most massive of these 
systems re ly instead upon artificially igniting a detona- 
tion front (P^^m or et al.|[2M0 2012 Raskin et al.||2014 


Moll et al. 2014 A propo sed model involving head- 


on collisions of white dwarfs (HawIe^etaL 2012 [Raskin 


et al.|200^|Kushnir et al.|2013[|Aznar-Siguan et al.|2013[ 

2014|) demonstrates robust detonations. However, it re- 
mams unclear whether head-on collisions, even in triple 
systems, match ob served SNe la rates and ®®Ni yields 
( Scalzo et al. [2014 1. 

h'or all prompt merger scenarios, it is crucial to recog¬ 
nize that the ®®Ni nucleosynthetic yield, and hence the 
peak luminosity, is closely correlated with the central 
densit y, or equivalently, th e mass of the white dwarf pri¬ 
mary ( Raskin et al.|[Ml4 ). Consequently, in the context 
of double-degenerate mergers, a deeper understanding of 
the most fundamental properties of SNe la, including 
both the absolute value and relatively small variance of 
their intrinsic luminosity, hinges crucially upon an expli¬ 
cation of their detonation mechanism. 


2. NUMERICAL SIMULATIONS 

We conducted a set of smoothed particle hydrodynam¬ 
ics (SPH) simulations of merging carbon-oxygen white 
dwarfs, with equal abundances of carbon and oxygen, 
and with masses of I.OMq and I.IMq. For such simula¬ 
tions we use d the SPH code employed in previous studies 
of this kind (Loren-Aguilar et al.|20lb), with a resolution 
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of 4 X 10® SPH particles. The overall time evolution of 
the merger process found in these simulations and the 
general appearance of the remnants are very similar to 
those found in other previous studies. 

At approximately the instant at which the maximum 
temperature is reached during the merger, the SPH con¬ 
ditions are remapped onto an Euler ian adaptive mesh 


using the FLASH4 code framework (Fryxell et al. 2000 
Dubey et al. 2009). For the model presented here, this 


remapping occurs 1.5 initial orbital periods, or 40 s, sub¬ 
sequent to tidal disruption. Throughout this paper, we 
refer to the time of remapping as t = 0 . 

FLASH solves the Euler equations of hydrodynamics 
using an unsplit hydrodynamics solver (Lee & Deane 
20091, and closed with an equation of state that includes 


contributions from blackbody radiation, ion s, and elec¬ 


trons of an arbitrary degree of degeneracy (Timmes & 


Swesty 2000). Nuclear burning is incorporated through 


the use ol a simplified 13 species alpha-chain network, 


which includes the effect of neutrino cooling (Timmes 
1999). The Poisson equation for self-gravity is solvecf 


both with a recently improved multipole solver (Couch 


et al. 2013), which treats gravitational potential at cell 


faces and centers the expansion at a square-density- 
weighted mean location, and also for comparis on pur¬ 


poses , with an oct-tree direct multigrid method (Ricker 


Like others (Pakmor et al.| 


20 T^ 

Raskin et al. 

2014 

Moll et al. 

2014), we do not 

tind that the initial 


at the onset of the three-dimensional FLASH simulation. 
However, unlike prior work, we choose not to initiate a 
detonation artihcially and instead follow the subsequent 
evolution of the system forward in time. Figure [l] dis¬ 
plays six snapshots of the overall evolution of the merger 
remnant. By t = 30 s, or roughly one initial orbital pe¬ 
riod, a pronounced m = 1 spiral mode develops, which 
persists for several dynamical times. At t = 108.98 s, a 
strong detonation driven by the spiral mode occurs self- 
consistently, giving rise to a supernova outburst. 

3. GRAVITATIONAL INSTABILITY AND SPIRAL MODES 

During the final inspiral of the binary white dwarf 
system, as the less massive secondary white dwarf ap¬ 
proaches the more massive primary, it is tidally dis¬ 
rupted, producing an accretion disk, with little mass 
ejected from the system. Conditions for thermonuclear 
detonation in degenerate carbon require hot t empera- 
tures T ~ 2 x 10^ K at densities p ~ 10^ g/cm® (Seiten- 


zahl et al. 2009a). In a sufficiently massive system, the 
primary has the requisite density for detonation, and the 
disk has the necessary temperature. However, for un¬ 
equal mass mergers, the tidal disruption of the secondary 
leaves the primary largely intact and unheated. Con¬ 
sequently, a key challenge faced by double-degenerate 
models is the need for some mechanism to raise the tem¬ 
perature of the cold primary {T ~ 10^ K) in order to 
produce a SNe la. If the mass accretion rate onto the 
central white dwarf is sufficiently rapid, it may mix hot, 
tidally disrupted disk material with the higher density 
white dwarf primary material, leading to unstable car¬ 
bon burning and a subsequent detonation, as pre vious 
studies of singl e-degenerate models have revealed (Jor- 
dan et aL|| 2012 ). 


Fundamental physical considerations determine the ba¬ 
sic physical state and evolution of the gas in the accretion 
disk, as a function of the primary mass Mi , the secondary 
mass M 2 , as well as the mass ratio q = M 2 /Mi. In par¬ 
ticular, the secondary mass is disrupted as its orbit ap¬ 
proaches the tidal radius i?tid = q~^^^R 2 , defined as the 
orbital separation distance at which the mean density of 
the primary, averaged over the entire spherical volume of 
the orbit, is equal to that of the secondary star, with an 
initial undeformed radius of i? 2 - 

We now show that the inner disk of a white dwarf 
merger is marginally gravitationally stable and subject 
to the development of non-axisymmetric spiral modes, 
based upon the physics of tidal disruption of the sec¬ 
ondary. We assume that the secondary is completely 
disrupted and extends radially from the tidal radius i?tid 
to large radii, while the primary remains intact. These 
assumptions are appropriate for mergers with unequal 
masses, which we verify a posteriori from numerical sim¬ 
ulation data. 

The sound speed Cg within the disk is then di¬ 
rectly determined from the thermal energy, Cs{r) = 
(/yGMi/r)®/^ in the non-degenerate ideal gas conditions 
prevalent in the disk, where 7 ~ 5/3 is the ratio of spe¬ 
cific heats, and / ^ 0.2 is the fraction of the maximum 
virial t emperature, as d etermined from numerical simu¬ 
lations (Zhu et al.|2013). In the inner disk, the rotational 
profile is nearly Keplerian, and the epicyclic frequency is 
K,{r) ~ H(r) ~ (GMi/r®)®/^. In addition, the surface 
gravity of the disk is determined by the surface density 
E(r) of the disk, which falls off as a power-law in ra¬ 
dius, 7 -“’*, with a relatively steep scaling of 77 cs 2.5 — 3.5, 
also d etermined from numerical simulations (Fryer et al. 
2010|). In the following, we presume a typical surface 

3 




density scaling F = (M2/27ri?/j^ 

These fundamental physical properties of the disk can 
then be combined to form the Toomre parameter Q = 
Csk/ttGE, which quantifies the degree to which the disk 
is subject to gravitational instability. Here, k is the 
epicyclic frequency, with k ^ H for a near-Keplerian 
disk. In particular, a classical analytic result yields that 
if Q < I, disks are subject to axisymmetric instability. 

Furthermore, models of massive accretion disks have 
demonstrated the existence of an eccentric gravitational 
instability even in disks that are otherwise gravitation¬ 
ally stable for disk masses exceeding a third that of the 
central stellar mass, provided that Q < 3 at corotation 
( Adams et al.|I989 Shu et al.|I990 ). The eccentric grav- 
itational instability gives rise to a characteristic m = 1 
spiral mode, which was originally studied in the context 
of massive star-forming accretion disks and has subse- 
quently been examined in a variety of other contexts 
( Hopkins fc Quataert||20I0 ). 

On the basis of these loasic physical considerations, 
Q{r) ~ q~^ir/Rtid) for binary white dwarf mergers, 
indicating marginal gravitational stability for a typical 
merger system. Furthermore, we note that this scal¬ 
ing of Q values is, remarkably, independent of the to¬ 
tal system mass. Consequently, we infer that many bi¬ 
nary dwarf mergers, including both sub-Chandraskehar 
and super-Chandrasekhar mass systems, with sufficiently 
large mass ratio q will generally be subject to spiral mode 
instabilities, including the eccentric gravitational insta- 
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bility. 

To elucidate the properties of the disk, we analyze the 
properties of the white dwarf merger system at t = 0 in 
figure The values and the power law scalings of the 
surface density S(r), sound speed Cs(r), angular veloc¬ 
ity fi(r), and Q(r) as determined from the full numerical 
simulation of the merger process are in very good agree¬ 
ment with our analytic estimates. In particular, Q ~ 1 in 
the inner accretion disk, confirming our prediction that 
the inner disk is marginally gravitationally stable. 

We further characterize the detonation mechanism in 
figure]^ This figure analyzes the structure of the n ierger 
at the instant o f detonation. As previous studies ( Blaes^ 
|& Hawley|1988|) have demonstrated, the spiral instability 
provides a mechanism that drives an inflow of hot mate¬ 
rial from the inner disk inward. The connection of the 
spiral to the funnel accretion flow is made clear in the 
bottom two panels of figure The hot matter from the 
tidally disrupted secondary is accreted inward through a 
funnel accretion flow generated by the spiral instability. 

4. DETONATION AND NUMERICAL TESTS 

The detonation is initiated in the degenerate material 
of the primary, where the density p = 6.7 x 10® g/cm®, 
and the temperature has been raised to T = 3.2 x 10® K 
just prior to runaway through the admixture of degener¬ 
ate white dwarf and hot disk material. Such a detona¬ 
tion arises due to the high virial temperature of the disk. 
Using a Salpeter IMF and the empirical initial -to-final 
white dwarf mass relation of|Kalirai et al. (2008), we es¬ 
timate that the fraction of WU binaries with a total mass 
greater than 2.1 Mq to be 2%. A much larger fraction of 
binaries will be susceptible to unstable helium burning 
(R. Kashyap et al. 2015, in preparation). 

The simulation undergoes a nuclear runaway, and the 
overpressure from the resulting nuclear energy injection 
initiates the detonation front. The minimum detonation 
length scale is 7 km, and is much smaller still for the 
carbon-burnin g zone, of cm order at t hese densities and 
temperatures (Seitenzahl et al.||2009b I. Both are clearly 
unresolved in all existing 3U simulations. However, we 
emphasize that the initiation is driven self-consistently 
by the spiral mode here without the need for an artificial 
detonation of any sort. Furthermore, an extensive set 
of verification tests demonstrates that the key proper¬ 
ties of the detonation are all robustly determined above 
a critical numerical resolution required to resolve spiral 
structure. A total of ^ 0.64 — O.72M0 of ®®Ni is synthe¬ 
sized during the detonation, corresponding to a normal¬ 
brightness SN la event. 

We have carried out extensive verification tests to con¬ 
firm the reliability and reproducibility of our simulations. 
In particular, we have carried out a convergence test, 
using higher and lower resolutions than our standard 
model. Additionally, because the use of the multipole 
gravity solver could be called into question in a strongly 
non-spherical flow, we repeated our simulations using an 
entirely independent multigrid gravity solver. The runs 
are summarized in table [l] where they are listed in order 
of decreasing resolution. Additionally, runs designated 
“mp” utilize the multipole solver, whereas “mg” runs 
employ multigrid. The domain size in all models is held 
fixed. 

At sufficiently high resolution, all models using multi- 


Table 1 

Runs in This Paper. Maximum finest resolution for each model is 
quoted along with the time of detonation (if any). Models that 
do not detonate list a dash (—) in place of a time of detonation, 
M(®®Ni), and max(E]jin). 


Run 

Resolution 

(km) 

^burn 

^det 

(S) 

M(56Ni) 

(Me) 

maxfEbin) 

(erg) 

1-mp 

68.3 km 

1 

89 s 

0.72 Me 

1.10 X 10®^ 

2-mpl 

136 km 

1 

112 s 

0.70 Me 

1.08 X 10®^ 

2-mp2 

136 km 

0.5 

112 s 

0.66 Me 

1.08 X 10®i 

2-mp3 

136 km 

0.3 

112 s 

0.64 Me 

1.08 X 10®i 

3-mp 

273 km 

1 

— 

— 

— 

4-mp 

546 km 

1 

— 

— 

— 

l-mg 

68.3 km 

1 

93.1 s 

0.72 Me 

1.06 X 10®^ 

2-mg 

136 km 

1 

107 s 

0.70 Me 

1.08 X 10®^ 

3-mg 

273 km 

1 

— 

— 

— 


grid and multipole produce nearly identical spiral struc¬ 
tures and detonations. More importantly, these mod¬ 
els demonstrate that the results are entirely insensitive 
to the choice of gravity solver at sufficiently high reso¬ 
lution. Further, these models demonstrate convergence 
at the level of ~ 3% in the quantity of ®®Ni produced. 
Lower resolution models (3-mp, 3-mg, 4-mp) fail to det¬ 
onate. Analysis of these models reveals that the spiral 
structure is much less pronounced in runs 3-mp and 3- 
mg and entirely absent from 4-mp. These findings are 
consistent with the interpretation that the detonation in 
these models is driven by the spiral structure generated 
in the accretion disk. 

Additionally, we have conducted a convergence study 
in the time domain by reducing the chosen timestep for 
a fixed spatial resolution. In particular, the timestep 
At” in the nth timestep is chosen to be the more re¬ 
strictive of the CFL timestep and the burn 

timestep e.g. Ar = min(At;;yd^„, At;)^^„). The 

burn timestep is in turn restricted to a value that pro¬ 
hibits the internal energy eint of a given cell from chang¬ 
ing by a fraction Cbum in a given timestep. That is, 
^^burn = At”"^Cburn(eint/enucAr"^), where e„uc is the 
specific rate of change of internal energy due to nuclear 
processes. 

Previous st udies of zero-impac t parameter colliding 
white dwarfs (Hawley et al. 2012) demonstrated a high 
degree of sensitivity to the burning time step criterion, 
through the parameter Cburm with ®®Ni yields varying 
by roughly 30%. These calculations suggested that the 
amount of ®®Ni generated through collisions is sensitive 
to the resolution of the simulation in the time domain. To 
check this dependence, we computed three sets of mod¬ 
els for our standard case, varying Cbum from I (effectively 
using the CFL condition only) to 0.5 and 0.3. We find 
a much stronger degree of sensitivity to the timestep se¬ 
lection than the spatial resolution, with the overall ®®Ni 
yield varying by slightly less than 10% over the range 
of models considered, with more highly-refined time do¬ 
main calculations producing less ®®Ni. However, over¬ 
all we find a smaller degree of sensitivity to the c hoice 
of the burning timestep than Hawley et al. (2012), and 
indeed there is a trend towards convergence evident in 
the more highly refined time domain runs. The use 
of the burning-adjusted timestep criterion dramatically 
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reduces the timestep relative to a calculation adopting 
the CFL condition only. In our simulations, the mini¬ 
mum burning-adjusted timesteps for Cbum = 0.5 and 0.3 
dropped to factors of 10“^ and 10“^, relative to the CFL 
timestep. Further reductions in Cbum to 0.2 dropped the 
timestep to a factor of 10“® relative to the CFL timestep, 
and proved to be computationally intractable. 

We note that the structure of the nuclear network in¬ 
volves a stiff reaction solver which evolves the nuclear re¬ 
actions at every At. However, the hydrodynamics is (by 
definition) not subcycled for an unsplit explicit solver. 
Consequently, the choice of the timestep directly deter¬ 
mines not so much the accuracy of the burning as the ac¬ 
curacy of hydrodynamical coupling to the nuclear burn¬ 
ing. In this capacity, it is not surprising that the head-on 
collision of white dwarfs, involving very strong shock- 
shock interactions on very small length scales, demon¬ 
strates a more sensitive coupling of the hydrodynamics 
to nuclear burning. 


5. CONCLUSION 

We have demonstrated that the innermost region of the 
accretion disk formed in the aftermath of a sufficiently 
massive white dwarf binary merger with a large mass 
ratio is accreted on the primary white dwarf due to spi¬ 
ral instabilities. The ensuing hot disk material accreted 
on to the central remnant drives a powerful thermonu¬ 
clear outburst, leading to a SN la explosion, without the 
need for artificial ignition. The spiral mechanism oper¬ 
ates on a dynamical timescale, about two to three orders 
of magnitude more rapidly than an alternative mecha¬ 
nism driven by the magnetorotational instability (Ji et al. 
2013|). Our results demonstrate the robustness of the 
double-degenerate scenario over a wider range of mass 
ratios than previously thought possible and shed light 
on the possible scenarios producing SN la, a topic of the 
major interest because of their application to cosmology. 
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Figure 1. Time evolution of the binary white dwarf system subsequent to merger. This figure presents the time evolution of the density 
in the mid-plane of the disk, shown on a logarithmic scale. The bottom two frames, taken at t = 108.98 and 109.88 s, are just prior to and 
just subsequent to the onset of detonation, respectively. The bottom two frames are zoomed-in by a factor of two in the x and y directions. 
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Figure 2. Properties of the accretion disk. Clockwise, from upper-left, this plot presents surface density S, sound speed Cs, angular 
velocity Q, and Toomre Q parameter versus cylindrical radius r as measured from the peak central density of the white dwarf merger, all 
on logarithmic-logarithmic axes, for the initial time t = 0 at the beginning of the simulations. The vertical dashed line indicates the radial 
location enclosing Mi = I.IMq, which demarcates the primary white dwarf on the left from the disk on the right. For comparison, we also 
indicate with a dashed line the location of r = -Rtid- The power law fits have indices of -3.25 for S, -0.96 for Cs, -1.60 for Q, and 0.92 for 
Toomre Q. The fits are in very good agreement with the simple scalings presented in the text. 
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Figure 3. Detailed view of the detonation. All four panels depict the first instant of the detonation tracked in the simulation, at f = 108.98 
s, corresponding to the bottom left panel of figure^ The upper-left panel shows the log(p), and shows the location of the detonation, near 
10^ g/cm^, as indicated by the white circle. The upper-right panel shows the log(T) where the peak temperature has risen to ~ 3 x 10® 
at the location of detonation. The bottom-left panel shows the divergence of the velocity field; positive values represent diverging flows 
and negative values indicate converging flows. The strongly negative values coincide with the spiral density waves of the first panel, and 
demonstrate the existence of shocks at the location of the spirals. The bottom-right panel shows the velocity streamlines in the disk, shaded 
by the radial velocity. The spiral structure directly gives rise to a funnel accretion inflow onto to the white dwarf primary, which mixes hot 
disk material into the merger and gives rise to the detonation. 
















